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We study the effect that the injection of a common source of noise has on the trajectories of chaotic 
systems, addressing some contradictory results present in the literature. We present particular 
examples of 1-d maps and the Lorenz system, both in the chaotic region, and give numerical evidence 
showing that the addition of a common noise to different trajectories, which start from different 
initial conditions, leads eventually to their perfect synchronization. When synchronization occurs, 
the largest Lyapunov exponent becomes negative. For a simple map we are able to show this 
phenomenon analytically. Finally, we analyze the structural stability of the phenomenon. 



The synchronization of chaotic systems has 
been the subject of intensive research in the last 
years. Besides its fundamental interest, the study 
of the synchronization of chaotic oscillators has a 
potential application in the field of chaos commu- 
nications. The main idea resides in the hiding of 
a message within a chaotic carrier generated by 
a suitable emitter. The encoded message can be 
extracted if an appropriate receiver, one which 
synchronizes to the emitter, is used. One of the 
conditions to be fulfilled in order to achieve syn- 
chronization is that the receiver and the emitter 
have very similar device parameters, hence mak- 
ing it very difficult to intercept the encoded mes- 
sage. Although the usual way of synchronizing 
two chaotic systems is by injecting part of the 
emitted signal into the receiver, the possibility 
of synchronization using a common random forc- 
ing has been also suggested. However, there have 
been some contradictory results in the literature 
on whether chaotic systems can indeed be syn- 
chronized using such a common source of noise 
and the issue has began to be clarified only very 
recently. In this paper we give explicit exam- 
ples of chaotic systems that become synchronized 
by the addition of Gaussian white noise of zero 
mean. We also analyze the structural stability of 
the phenomenon, namely, the robustness of the 
synchronization against a small mismatch in the 
parameters of the chaotic sender and receiver. 



I. INTRODUCTION 

One of the most surprising results of the last decades 
in the field of stochastic processes has been the discover- 
ing that fluctuation terms (loosely called noise) can ac- 
tually induce some degree of order in a large variety of 
non-linear systems. The first example of such an effect 
is that of stochastic resonance by which a bistable 
system responds better to an external signal (not nec- 



essarily periodic) under the presence of fluctuations, ei- 
ther in the intrinsic dynamics or in the external input. 
This phenomenon has been shown to be relevant for some 
physical and biological systems described by nonlinear 
dynamical equations |^-^|. Other examples in purely 
temporal dynamical systems include phenomena such as 
noise-induced transition s |H ] , noise-induced transport , 
coherence resonance &-ITW, etc. In extended systems, 
noise is known to induce a large variety or ordering ef- 
fects [p"2| , such as pattern formation 1^,0 , phase tran- 
sitions~|l5|-|l8|l , phase separation flt^ , [20|| , spatiotemporal 



noise-sustained structures 
doubly stochastic resonance |25|1, amongst many 



stochastic resonance 

m 

others. All these examples have in common that some 
sort of order appears only in the presence of the right 
amount of noise. 

There has been also some recent interest on the in- 
terplay between chaotic and random dynamics. Some 
counterintuitive effects such as coherence resonance, or 
the appearance of a quasi-periodic behavior, in a chaotic 
system in the presence of noise, have been found recently 
P6|| . The role of noise in standard synchronization of 
chaotic systems has been considered in ^7j,^8|, as well 
as the role of noise in synchronizing non-chaotic systems 
p9| , pp[ . In this paper we address the different issue of 
synchronization of chaotic systems by a common random 
noise source, a topic that has attracted much attention 
recently. The accepted result is that, for some chaotic 
systems, the introduction of the same noise in indepen- 
dent copies of the systems could lead (for large enough 
noise intensity) to a common collapse onto the same tra- 
jectory, independently of the initial condition assigned to 
each of the copies. This synchronization of chaotic sys- 
tems by the addition of random terms is a remarkable and 
counterintuitive effect of noise and although some clarify- 
ing papers have appeared recently, still some contradic- 
tory results exist for the existence of this phenomenon 
of noise-induced synchronization. It is the purpose of 
this paper to give further analytical and numerical evi- 
dence that chaotic systems can synchronize under such 
circumstances and to analyze the structural stability of 
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the phenomenon. Moreover, the results presented here 
clarify the issue, thus opening directions to obtain such 
a synchronization in electronic circuits, for example for 
encryption purposes. Common random noise codes have 
been used in spread spectrum communication since a long 
time ago f3f|] . The main idea is to mix a information data 
within a noisy code. At the receiver, the information is 
recovered using a synchronized replica of the noise code. 
More recently, the use of common noise source has been 
also proposed as a useful technique to improve the en- 
cryption of a key in a communication channel p2| . 

The issue of ordering effect of noise in chaotic systems 
was considered already at the beginning of the 80's by 
Matsumoto and Tsuda Q who concluded that the in- 
troduction of noise could actuall y m ake a system less 
chaotic. Later, Yu, Ott and Chen |34| studied the transi- 
tion from chaos to non-chaos induced by noise. Synchro- 
nization induced by noise was considered by Fahy and 
Hamman |35| who showed that particles in an external 
potential, when driven by the same random forces, tend 
to collapse onto the same trajectory, a behavior inter- 
preted as a transition from chaotic to non-chaotic behav- 
iors. The same system has been studied numerically and 
analytically []36|-|38| . Pikovsky |39| analyzed the statistics 
of deviations from this noise-induced synchronization. A 
paper that generated a lot of controversy was that of 
Maritan and Banavar jfo). These authors analyzed the 
logistic map in the presence of noise: 



x n +i = 4a;„(l - x n ) + £n 



(1) 



where £ n is the noise term, considered to be uniformly 
distributed in a symmetric interval [— W, +W]. They 
showed that, if W was large enough (i.e. for a large noise 
intensity) two different trajectories which started with 
different initial conditions but used otherwise the same 
sequence of random numbers, would eventually coincide 
into the same trajectory. The authors showed a similar 
result for the Lorenz system (see section III). This result 
was heavily criticized by Pikovsky ]4l| ] who proved that 
two systems can synchronize only if the largest Lyapunov 
exponent is negative. He then argued that the largest 
Lyapunov exponent of the logistic map in the presence 
of noise is always positive and concluded that the syn- 
chronization was, in fact, a numerical effect of lack of 
precision of the calculation. The analysis of Pikovsky 
was confirmed by Longa et al. jf|| who studied the logis- 
tic map with arbitrary numerical precision. The criterion 
of negative Lyapunov exponent has also been shown to 
hold for other types of synchronization of chaotic sys- 
tems and Zhou and Lai |43| ] noticed that previous results 
by Shuai, Wong and Cheng [Q showing synchronization 
with a positive Lyapunov exponent were again an artifact 
of the limited precision of the calculation. 

In addition to the above criticisms, Herzel and Fre- 
und |^5| and Malescio [[l6j pointed out that the noise 
used to simulate Eq.(Q) and the Lorenz system in ||(]] 



is not really symmetric. While the noise in the Lorenz 
system is non-symmetric by construction, in the case of 
the map, the non-zero mean arises because the require- 
ment x n £ (0, 1), Vn, actually leads to discard the values 
of the random number which would induce a viola- 
tion of such condition. The average value of the random 
numbers which have been accepted is different from zero, 
hence producing an effective biased noise, i.e. one which 
does not have zero mean. The introduction of a non-zero 
mean noise means that the authors of ^] were altering 
essentially the properties of the deterministic map. Fur- 
thermore, Cade and Bassu jf7| argued that the synchro- 
nization observed by Maritan and Banavar is due to the 
fact that the bias of the noise leads the system to a non- 
chaotic fixed point. With only this basis, they concluded 
that a zero-mean noise can never lead to synchronization 
in the Lorenz system. The same conclusion was reached 
by Sanchez et al. (48| who studied experimentally a Chua 
circuit and concluded that synchronization by noise only 
occurs if the noise does not have a zero mean. The same 
conclusion is obtained in flf9| by studying numerically a 
single and an array of Lorenz models, and in ]50[ | from 
experiments in an array of Chua circuits with multiplica- 
tive colored noise. Therefore, from these last works, a 
widespread belief has emerged according to which it is 
not possible to synchronize two chaotic systems by in- 
jecting the same noisy unbiased, zero-mean, signal to 
both of them. 

Contrary to these last results (but in agreement with 
the previously mentioned results ]32]-|39[|), Lai and Zhou 
f3lf have shown that some chaotic maps can indeed be- 
come synchronized by additive zero-mean noise. A sim- 
ilar result has been obtained by Loreto et al. f5^| , and 
by Minai and Anand |32|,^3 54 1, in the case where the 
noise appears parametrically in the map. The implica- 
tions to secure digital communications have been consid- 
ered in IHJHUm , and an application to ecological dy- 
namics in fluid flows is presented in fl56| . An equivalent 
result about the synchronization of Lorenz systems using 
a common additive noise has been shown by the authors 
of the present paper in jSTy . The actual mechanism that 
leads to synchronization has been explained by Lai and 
Zhou [^lj, see also |5q] . As Pikovsky Q required, syn- 
chronization can only be achieved if the Lyapunov expo- 
nent is negative. The presence of noise allows the system 
to spend more time in the "convergence region" where 
the local Lyapunov exponent is negative, hence yielding 
a global negative Lyapunov exponent. This argument 
will be developed in more detail in section [H], where an 
explicit calculation in a simple map will confirm the anal- 
ysis. The results of Lai and Zhou have been extended to 
the case of coupled map lattices (5^j where Pikovsky's 
criterion has been extended for spatially extended sys- 
tems. 

In this paper we give further evidence that it is pos- 
sible to synchronize two chaotic systems by the addition 
of a common noise which is Gaussianly distributed and 
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not biased. We analyze specifically some 1-d maps and 
the Lorenz system, all in the chaotic region. The nec- 
essary criterion introduced in Ref. |yj and the general 
arguments of |5l) are fully confirmed and some heuris- 
tic arguments are given about the general validity of our 
results. 

The organization of the paper is as follows. In sec- 
tion H we present numerical and analytical results for 
some 1-d maps, while sectianJjTJ studies numerically the 
Lorenz system. In section |IV| we analyze the structural 
stability of the phenomenon, i.e. the dependence of the 
synchronization time on the parameter mismatch. Fi- 
nally, in section [V| we present the conclusions as well as 
some open questions relating the general validity of our 
results. 



of both variables becoming themselves identical to the 
noise term. 



II. RESULTS ON MAPS 

The first example is that of the map: 

X n +1 = F(x n ) = f[x n ) + e£n 



(2) 



where is a set of uncorrelated Gaussian variables of 
zero mean and variance 1. As an example, we use explic- 
itly 



exp 



a; -0.5 



(3) 



Studying the convergence or divergence of trajectories of 
Eq. (^) starting from different initial conditions under 
the same noise £„ is equivalent to analyzing the converge 
or divergence of trajectories from two identical systems 
of the form (^|) driven by the same noise. We plot in 
Fig.(l) the bifurcation diagram of this map in the noise- 
less case. We can see the typical windows in which the 
system behaves chaotically. The associated Lyapunov ex- 
ponent, A, is positive in these regions. For instance, for 
lu = 0.3 (the case we will be considering throughout the 
paper) it is A » 0.53. In Fig. (2) we observe that the 
Lyapunov exponent becomes negative for most values of 
lu for large enough noise level e. Again for lu = 0.3 and 
now for e = 0.2 it is A = —0.17. A positive Lyapunov 
exponent in the noiseless case implies that trajectories 
starting with different initial conditions, but using the 
same sequence of random numbers {£«}, remain different 
for all the iteration steps. In this case, the corresponding 
synchronization diagram shows a spread distribution of 
points (see Fig. (3a)). However, when moderate levels of 
noise (e > 0.2) are used, A becomes negative and trajec- 
tories starting with different initial conditions, but using 
the same sequence of random numbers, synchronize per- 
fectly, see the synchronization diagram in Fig. (3b). Ob- 
viously, the noise intensity in the cases shown is not large 
enough such as to be able to neglect completely the de- 
terministic part of the map. Therefore, the synchroniza- 
tion observed does not trivially appear as a consequence 




FIG. 1. Bifurcation diagram of the map given by Eqs. 
and (^) in the absence of noise terms. 
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FIG. 2. Lyapunov exponent for the noiseless map (e = 0, 
continuous line) and the map with a noise intensity e = 0.1 
(dotted line) and e = 0.2 (dot-dashed line). 
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FIG. 3. Plot of two realizations of of the map 

given by Eqs. (^|) and (^) with to = 0.3. Each realization 
consists of 10,000 points which have been obtained by itera- 
tion of the map starting in each case from a different initial 
condition (100,000 initial iterations have been discarded and 
are not shown). In figure (a) there is no noise, e = and 
the trajectories are independent of each other. In figure (b) 
we have used a level of noise e = 0.2, producing a perfect 
synchronization (after discarding some initial iterations). 

According to jyj, convergence of trajectories to the 
same one, or lack of sensitivity to the initial condition, 
can be stated as negativity of the Lyapunov exponent. 
The Lyapunov exponent of the map (^) is defined as 

1 N 

N— >oo /V z — ' 

i=l 

It is the average of (the logarithm of the absolute value 
of) the successive slopes F' found by the trajectory. 
Slopes in [—1,1] contribute to A with negative values, in- 
dicating trajectory convergence. Larger or smaller slopes 
contribute with positive values, indicating trajectory di- 
vergence. Since the deterministic and noisy maps satisfy 
F' = f one is tempted to conclude that the Lyapunov ex- 
ponent is not modified by the presence of noise. However, 
there is noise-dependence through the trajectory values 
Xj, i = 1, N. In the absence of noise, A is positive, in- 
dicating trajectory separation. When synchronization is 
observed, the Lyapunov exponent is negative, as required 
by the argument in J4l| . 

Notice that this definition of the Lyapunov exponent 
assumes a fixed realization of the noise terms, and it is 
the relevant one to study the synchronization phenomena 
addressed in this paper. One could use alternative def- 
initions ]52| . For instance, if one considers the coupled 
system of both the x variable and the noise generator 
producing £, then the largest Lyapunov exponent of the 



composed system is indeed positive (and very large for 
a good random number generator). This simply tells us 
that there is a large sensitivity to the initial condition 
of the composed system (x, £) as shown by the fact that 
a change of the seed of the random number generator 
completely changes the sequence of values of both £ and 
x. We consider in this paper the situation described by 
definition (Q) with fixed noise realization. 

By using the definition of the invariant measure on the 
attractor, or stationary probability distribution P st (x), 
the Lyapunov exponent can be calculated also as 

\={log\F'(x)\) = (log|/'(x)|) = J P st (x)log\f'(x)\dx 

(5) 

Here we see clearly the two contributions to the Lyapunov 
exponent: although the derivative f'(x) does not change 
when including noise in the trajectory, the stationary 
probability does change (see Fig. 4), thus producing the 
observed change in the Lyapunov exponents. Synchro- 
nization, then, can be a general feature in maps, such 
as (|^), which have a large region in which the deriva- 
tive |/'(x)| is smaller than one. Noise will be able to 
explore that region and yield, on the average, a negative 
Lyapunov exponent. This is, basically, the argument de- 
veloped in |nj. 



(a) 




FIG. 4. Plot of the stationary distribution for the map 
given by Eqs.(^) and (^) with lu = 0.3 in the (a) deterministic 
case e = 0, and (b) the case with noise along the trajectory, 
e = 0.2. 

In order to make some analytical calculation that can 
obtain in a rigorous way the transition from a positive 
to negative Lyapunov exponent, let us consider the map 
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given by Eq. (||) and 



/(*) = 



a(l 


- exp(l + 


x)) if x < -1 


-2 


- 2x 


if x £ (-l,--5) 


2x 




if x G (-.5, .5) 


2- 


2x 


if x € (.5,1) 


a(- 


1 + exp(l 


- x)) if a; > 1 



(0) 



with < a < 1. This particular map, based in the tent 
map |^3(|, has been chosen just for convenience. The fol- 
lowing arguments would apply to any other map that in 
the absence of noise takes most frequently values in the 
region with the highest slopes, but which visits regions of 
smaller slope when noise is introduced. This is the case, 
for example, of the map (3). In the case of (||), the val- 
ues given by the deterministic part of the map, after one 
iteration from arbitrary initial conditions, fall always in 
the interval (—1, 1). This is the region with the highest 
slope \F'\ = 2. In the presence of noise the map can take 
values outside this interval and, since the slopes encoun- 
tered are smaller, the Lyapunov exponent can only be 
reduced from the deterministic value. To formally sub- 
stantiate this point, it is enough to recall the definition of 
Lyapunov exponent Q): an upper bound for |i*"(a;)| is 2, 
so that a bound for A is immediately obtained: A < In 2. 
Equality is obtained for zero noise. 

The interesting point about the map (||) and similar 
ones is that one can demonstrate analytically that A can 
be made negative. The intuitive idea is that it is enough 
to decrease a in order to give arbitrarily small values to 
the slopes encountered outside (—1, 1), a region accessible 
only thanks to noise. To begin with, let us note that 
\F'(x)\ = 2 if x e (-1, 1), and \F'(x)\ < a if |x| > 1, so 
that an upper bound to (0) can be written as 



A < lim 

N^oc 



^M2. 
N 



—jQ In a ) = /»/ In 2 + i>,,]]]<i 



In 2 -p Ln(2/a) 



(7) 



Ni/N and No /N are the proportion of values of the map 
inside I = (0, 1) and outside this interval, respectively, 
and we have used that as N — > oo they converge to pi 
and po, the invariant measure associated to I and to 
the rest of the real line, respectively (pi + po = 1)- A 
sufficient condition for x n +\ = f(x n ) + e£ n to fall outside 
I is that |£„| > 2/e. Thus, po — Probability(|x„ + i | > 
1) > Probability( |£„ | > 2/e) = crfc(V2/e) = T, where 
we have used the Gaussian character of the noise. In 
consequence, one finds from (^) 



A < ln2-Tln(2/a). 



(8) 



The important point is that T — erfc(v2/e) is indepen- 
dent on the map parameters, in particular on a. Thus, 
(||) implies that by decreasing a the value of A can 
be made as low as desired. By increasing e such that 
T > ln2/ln(2/a), A will be certainly negative. Thus we 
have shown analytically that strong enough noise will al- 
ways make negative the Lyapunov exponent of the map 



(||) and, accordingly, it will induce yield "noise-induced 
synchronization" in that map. 



III. THE LORENZ SYSTEM 

In this section we give yet another example of noise 
induced synchronization. We consider the well known 
Lorenz |6j| model with additional random terms of the 
form (|C[: 



±=p(y- x) 

y = —xz + rx — y - 

z — xy — bz 



(9) 



£ is white noise: a Gaussian random process of mean zero, 
(£(*)) = and delta-correlated, (£(*)£(*')) = S(t - t'). 
We have used p — 10, b — 8/3 and r = 28 which, in 
the deterministic case, e = are known to lead to a 
chaotic behavior (the largest Lyapunov exponent is A w 
0.9 > 0). As stated in the introduction, previous results 
seem to imply that synchronization is only observed for a 
noise with a non-zero mean. However, our results show 
otherwise. 

We have integrated numerically the above equations 
using the stochastic Euler method |6^]. Specifically, the 
evolution algorithm reads: 



x(t + At) = x(t) + At [p(y(t) - x(t))] 
y(t + At) = y(t) + At \-x{t)z{t) + rx{t) 

+ eVAtg(t) 
z(t + At) = z(t) + At [x{t)y(t) - bz(t)] 



y(t)} (io) 



The values of g(t) are drawn at each time step from an 
independent Gaussian distribution of zero mean and vari- 
ance one and they have been generated by a particularly 
efficient algorithm using a numerical inversion technique 
f33f . The time step used is At = 0.001 and simulations 
range typically for a total time of the order of t = 10 4 
(in the dimcnsionlcss units of the Lorenz system of equa- 
tions). The largest Lyapunov exponent has been com- 
puted using a simultaneous integration of the linearized 
equations p4| . For the deterministic case, trajectories 
starting with different initial conditions are completely 
uncorrelated, see Fig. (5a). This is also the situation for 
small values of e. However, when using a noise intensity 
e = 40 the noise is strong enough to induce synchroniza- 
tion of the trajectories. Again, the presence of the noise 
terms forces the largest Lyapunov exponent to become 
negative (for e = 40 it is A w —0.2). As in the exam- 
ples of the maps, after some transient time, two differ- 
ent evolutions which have started in completely different 
initial conditions synchronize towards the same value of 
the three variables (see Fig. (5b) for the z coordinate). 
Therefore, these results prove that synchronization by 



5 



common noise in the chaotic Lorenz system does occur 
for sufficiently large noise intensity. This result contra- 
dicts previous ones in the literature jh|[|8| . The main dif- 
ference with these papers is in the intensity of the noise: 
it has to be taken sufficiently large, as here, in order to 
observe synchronization. Notice that although the noise 
intensity is large, the basic structure of the "butterfly" 
Lorenz attractor remains present as shown in Fig. (6). 
Again, this result shows that, although the noise inten- 
sity used could be considered large, the synchronization 
is rather different from what would be obtained from a 
trivial common synchronization of both systems to the 
noise variable by neglecting the deterministic terms. 



IV. STRUCTURAL STABILITY 



An important issue concerns the structural stability of 
this phenomenon, in particular how robust is noise syn- 
chronization to small differences between the two systems 
one is trying to synchronize. Whether or not the syn- 
chronization of two trajectories of the same noisy Lorenz 
system (or of any other chaotic system) observed here, 
equivalent to the synchronization of two identical sys- 
tems driven by a common noise, could be observed in 
the laboratory, depends on whether the phenomenon is 
robust when allowing the two Lorenz systems to be not 
exactly equal (as they can not be in a real experiment). 
If one wants to use this kind of stochastic synchroniza- 
tion in electronic emitters and receivers (for instance, as a 
means of encryption) one should be able to determine the 
allowed discrepancy between circuits before the lack of 
synchronization becomes unacceptable. Additional dis- 
cussions on this issue may be found in p9| , p2| , |65[ . 
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FIG. 5. Plot of two realizations z' 1 ', of the Lorenz 
system Eqs.(|) with p = 10, b = 8/3 and r = 28. Each 
plotted realization starts from a different initial condition and 
consists of an initial warming up time of t — 12000 (not shown 
in the figure) and runs for a time t = 600 in the dimensionless 
units of the Lorenz system of equations. Panel (a) shows the 
deterministic case (e = 0) and panel (b) shows the results for 
e = 40. Notice the perfect synchronization in case (b). 
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FIG. 6. "Butterfly" attractor of the Lorenz system in the 
cases (a) of no noise e = 0, and (b) e — 40 using the same 
time series as in figure 5. 

We consider the following two maps forced by the same 
noise: 

x n +i = f(x„) +£„ (11) 
Vn+i = g{y n ) + £n (12) 

Linearizing in the trajectory difference u n — y n — x n , 
assumed to be small, we obtain 

u n +i = g'(x n )u n + g(x n ) - f(x n ) ee g'{x n )u n + A(x n ) 

(13) 

We have defined A(x) ee g{x) — f(x), and we are inter- 
ested in the situation in which the two systems are just 
slightly different, for example, because of a small param- 
eter mismatch, so that A will be small in some sense 
specified below. 

Iteration of (|ll]) leads to the formal solution: 

n— 1 

Un = M(n - l,0)uo + ^2 M( n - l,m + l)A(x m ) (14) 

m— 

We have defined M(j, i) — n*=i g'(x k ), and M[i — 1, i) = 
1. An upper bound on ([lj) can be obtained: 



G 



n-1 



Kl 2 < \M(n - 1,0)\ 2 \u \ 2 +J2 |M(n-l,m + l)| 2 |A(a 

(15) 



m=0 



The first term in the r.h.s. is what would be obtained 
for identical dynamical systems. We know that M(n — 
1,0)—* e Xn as n — > oo, where A is the largest Lyapunov 
exponent associated to (|l^). We are interested in the 
situation in which A < 0, for which this term vanishes 
at long times. Further analysis is done first for the case 
in which A(x) is a bounded function (or a; is a bounded 
trajectory with A continuous). In this situation, there is 
a real number /i such that |A(x m )| < /i. We then get: 



< /i 2 \M(n- l,m + l) 



(16) 



m=0 



an unequality valid for large n. Let us now define K = 
maxj |s(a;)|, the maximum slope of the function g{x). A 
trivial bound is now obtained as: 



\u„\ 2 < [i 2 



2)1 



1-K 
1-K 2 



(17) 



This can be further improved in the case K < 1, where 
we can write: 



\u n \ 2 < /i 2 



1 



1-K 2 



(18) 



As a consequence, differences in the trajectories remain 
bounded at all iteration steps n. Since, according to the 
definition (||), In if is also an upper bound for the Lya- 
punov exponent for all values of e and, in particular, for 
the noiseless map, e = 0, this simply tells us that if the 
deterministic map is non-chaotic, then the addition of a 
common noise to two imperfect but close replicas of the 
map will still keep the trajectory difference within well 
defined bounds. The situation of interest here, however, 
concerns the case in which a negative Lyapunov exponent 
arises only as the influence of a sufficiently large noise 
term, i.e. the deterministic map is chaotic and K > 1. 

lq) contains products of 



er than 1. It is still true 



In this case, the sum in Eq. 
slopes which are larger or small 
that the terms in the sum for large value of n — m can 
be approximated by M(n — 1, m + 1) w e (,n-m—i)X anc ^ 
considering this relation to be valid for all values of n, m, 
we would get: 



|Af(n-l,m+l)| 2 «e 2A ™ ^ 



-2A(m+l) 



1-e 



2An 



m=0 



m=0 



l_ e 2A 

(19) 



and, thus, at large n: 

KI 2 <M 2 (i- e 2A )^ 



(20) 



It can happen, however, that the product defining M (n — 
1, m + 1) contains a large sequence of large slopes g'(xi). 



TJiese terms (statistically rare) will make the values of 
l }li„| to violate the above bound at sporadic times. Analy- 
sis of the statistics of deviations from synchronization was 
carried out in |3j| . Although for A < the most probable 
deviation is close to zero, power-law distributions with 
long tails are found, and indeed its characteristics are 
determined by the distribution of slopes encountered by 
the system during finite amounts of time, or finite-time 
Lyapunov exponents, as the arguments above suggest. 
Therefore, we expect a dynamics dominated by relatively 
large periods of time during which the difference between 
trajectories remains bounded by a small quantity, but in- 
termittently interrupted by bursts of large excursions of 
the difference. This is indeed observed in the numerical 
simulations of the maps defined above. This general pic- 
ture is still valid even if |A(x)| is not explicitly bounded. 

We have performed a more quantitative study for the 
case in which two noisy Lorenz systems with different 
sets of parameters, namely: 



x\ 



■Pi{yi - xx) 

■■ -xizi + r x xi -yi + e£ 
- xiyi - b\%\ 



(21) 



and 



X2 

m 

Z2 



P2(V2 - X 2 ) 

-x 2 z 2 + r 2 x 2 - V2 + e£ 
X2V2 - b 2 z 2 



(22) 



are forced by the same noise In order to discern 

the effect of each parameter separately, we have varied 
independently each one of the three parameters, (p, 6, r), 
while keeping constant the other two. The results are 
plotted in Fig. 7. In this figure we plot the percentage of 
time in which the two Lorenz systems are still synchro- 
nized with a tolerance of 10%. This means that trajec- 
tories are considered synchronized if the relative differ- 
ence in the z variable is less than 10%. According to the 
general discussion for maps, we expect departures from 
approximate synchronization from time to time. They 
are in fact observed, but from Fig. 7 we conclude that 
small variations (of the order of 1%) still yield a syn- 
chronization time of more than 85%. In Fig. 8 we show 
that the loss of synchronization between the two systems 
appears in the form of bursts of spikes whose amplitude 
is only limited by the size of the attractor in the phase 
space. Moreover, it can be clearly seen in the same figure 
that large (but infrequent) spike amplitudes appear for 
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arbitrarily small mismatch. 




FIG. 7. Percentage of time that two slightly dissimilar 
Lorenz systems, subjected to the same noise, remain synchro- 
nized (up to a 10% discrepancy in the z variable). For one of 
the two systems we fix pi = 10, bi = 8/3 and n = 28 while 
for the other we vary systematically one of the parameters 
keeping the other two constant: in panel (a) the parameter 
62 varies , in panel (b) the parameter ri varies and in panel 
(c) the parameter P2 varies. Notice that the percentage of 
synchronization time is still higher than 85% if the relative 
difference between the parameters is less than 1%. 



In the realm of synchronization of chaotic oscillators, 
two different types of analogous intermittent behavior 
have been associated also to the fluctuating character of 
the finite-time conditional Lyapunov exponents as above. 
One is on-off intermitency p6| where the synchroniza- 
tion manifold is sligthly unstable on average but the fi- 
nite time Lyapunov exponent is negative during relatively 
long periods of time. In the other one, named bubbling 
|?7fl , the synchronization is stable on average but the lo- 
cal conditional Lyapunov exponent becomes occasionally 
positive. While in the former case bursting always occurs 
due to the necessarily imperfect initial synchronization, 
in the latter it is strictly a consequence of the mismatch 
of the entraining systems. In this sense, the behavior 
reported in the preceding paragraph should be consid- 
ered as a manifestation of bubbling in synchronization 
by common noise. 



40 
20 
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FIG. 8. Time evolution of the difference between two 
trajectories Zi and 22 corresponding to two Lorenz systems 
driven by the same noise but with a small mismatch in the r 
parameter: panel (a) n = 28 and 7-2 = 0.99999 x n and in 
panel (b) r-z — 0.999 x n. Notice that although the synchro- 
nization time diminishes with increasing parameter mismatch, 
the maximum absolute difference between the two variables 
does not depend on the mismatch. 



V. CONCLUSIONS AND OPEN QUESTIONS 

In this paper we have addressed the issue of synchro- 
nization of chaotic systems by the addition of common 
random noises. We have considered three explicit ex- 
amples: two 1-d maps and the Lorenz system under the 
addition of zero-mean, Gaussian, white noise. While the 
map examples confirm previous results in similar maps, 
and we have obtained with them analytical confirmation 
of the phenomenon, the synchronization observed in the 
Lorenz system contradicts some previous results in the 
literature. The reason is that previous works considered 
noise intensities smaller than the ones we found neces- 
sary for noise-synchronization in this system. Finally, 
we have analyzed the structural stability of the observed 
synchronization. In the Lorenz system, synchronization 
times larger than 85% (within an accuracy of 10%) can 
still be achieved if the parameters of the system are al- 
lowed to change in less than 1%. 

It is important to point out that noise-induced syn- 
chronization between identical systems subjected to a 
common noise is equivalent to noise induced order, in 
the sense that the Lyapunov exponent defined in (4) be- 
comes negative in a single system subjected to noise. One 
can ask whether the state with negative Lyapunov expo- 
nent induced by noise may be still be called 'chaotic' or 
not. This is just a matter of definition: if one defines 
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chaos as exponential sensibility to initial conditions, and 
one considers this for a fixed noise realization, then the 
definition of Lyapunov exponent implies that trajecto- 
ries are not longer chaotic in this sense. But one can 
also consider the extended dynamical system containing 
the forced one and the noise generator (for example, in 
numerical computations, it would be the computer ran- 
dom number generator algorithm). For this extended 
system there is strong sensibility to initial conditions in 
the sense that small differences in noise generator seed 
leads to exponential divergence of trajectories. In fact, 
this divergence is at a rate given by the Lyapunov ex- 
ponent of the noise generator, which approaches infin- 
ity for a true Gaussian white process. Trajectories in 
the noise-synchronized state are in fact more irregular 
than in the absence of noise, and attempts to calculate 
the Lyapunov exponent just from the observation of the 
time series will lead to a positive and very large value, 
since it is the extended dynamical system the one which 
is observed when analyzing the time series |6^] (typically 
such attempts will fail because the high dimensionality 
of good noise generators, ideally infinity, would put them 
out of the reach of standard algorithms for Lyapunov ex- 
ponent calculations). Again, whether or not to call such 
irregular trajectories with just partial sensibility to initial 
conditions 'chaotic' is just a matter of definition. More 
detailed discussion along these lines can be found in [B2| . 

There remain still many open questions in this field. 
They involve the development of a general theory prob- 
ably based in the invariant measure, that could give us 
a general criterion to determine the range of parame- 
ters (including noise levels) for which the Lyapunov ex- 
ponent becomes negative, thus allowing synchronization. 
In this work and similar ones, the word synchronization 
is used in a very restricted sense, namely: the coinci- 
dence of asymptotic trajectories. This contrasts with the 
case of interacting periodic oscillations where a more gen- 
eral theory of synchronization exists to explain the phe- 
nomenon of non trivial phase locking between oscillators 
that individually display very different dynamics. Indica- 
tions of the existence of analogue non trivial phase lock- 
ing have been reported for chaotic attractors |59| . There 
a "phase" with a chaotic trajectory defined in terms of a 
Hilbert transform is shown to be synchronizable by ex- 
ternal perturbations in a similar way as it happens with 
periodic oscillators. Whether or not this kind of general- 
ized synchronization can be induced by noise is, however, 
a completely open question. Last, but not least, it would 
be also interesting to explore whether analogs of the re- 
cently reported synchronization of spatio-temporal chaos 
|70| , |n| may be induced by noise. 
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